AD-A114  472  WISCONSIN  UNIV-MADISON  MATHEMATICS  RESEARCH  CENTER  F/6 


A  NEW  APPROACH  TO  EULER  SPLINES  II.(U) 

DEC  81  I  J  SCHOENBERG  DAAG29-80-C-0041 


END 


6 'l  £2 


OTIC  FILE  COPY 


B 


MRC  Technical  Summary  Report  #  2313 


A  NEW  APPROACH  TO  EULER  SPLINES  II 


I.  J.  Schoenberg 


Mathematics  Research  Center 
University  of  Wisconsin— Madison 
610  Walnut  Street 
Madison,  Wisconsin  53706 


December  1981 


(Received  October  28,  1981) 


Sponsored  by 

U.  S.  Army  Research  Office 
P.  O.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


Approved  for  public  release 
Distribution  unlimitod 


s 


DTK 

,  ELECTS 
MAY  1 8  888 


8Z  Ob 


Oa  8 


DT4C 


J 


UNIVERSITY  OF  WISCONSIN  -  MADISON 
MATHEMATICS  RESEARCH  CENTER 

A  NEW  APPROACH  TO  EULER  SPLINES  II 

I.  J*  Schoenberg 

Technical  Summary  Report  #2313 

December  1981 

ABSTRACT 


Accession  For 


NTIS  GRAAI 
DTIC  TAB 
Unannounced 
Justification- 


¥ 

□ 


By- 


Distribution/ 


Availability  Codes 
[Avail  and/or 
Special 


Let  t  G  C  be  a  constant  which  is  not  real  and  negative,  t  +  0,  t  +  1. 

In  [4]  it  was  shown  that  the  exponential  Euler  splines  Sn(x)  *  Sn(xjt),  *s 

introduced  by  the  author  in  [2],  can  be  obtained  by  the  following  recursive 

procedure:  We  define  S^x)  as  the  cardinal  linear  spline  interpolating  the 

v 

biinfinite  sequence  (t  )  (v  G  2)  and  then  defined  recursively 

x.i  1 

Sn(x)  =  /  1  Sn  1  (u)du//2  Sn  (u)du  (n  -  2,3,...)  . 

X~  2  "2 


In  [4]  this  was  used  to  derive  all  the  known  properties  of  these  splines  and 
also  some  new  ones.  In  the  present  short  note,  written  for  Euler's 
bicentenial  of  1983,  we  just  show  that  the  resulting  Sn(x)  are  identical 
with  the  splines  defined  in  [2]. 
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SIGNIFICANCE  AND  EXPLANATION 


Inthe  first paper  [4]  on  the  subject  the  new  recursive  approach  to  the 
exponential  Euler  splines  was  used  to  derive  the  known  properties^  and  some 

new  ones,  of  these  most  attractive  among  cardinal  splines.  In  the  present 

) 

short  note,  Witten  for  Euler's  bicentenial  of  1983,  we  }net  show  that  the  new 
recursive  construction  leads  precisely  to  the  exponential  Euler  splines*  as  w 
^Introduce*!  in  [2)«  v 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report* 


A  NEW  APPROACH  TO  EULER  SPLINES  II 


I.  J.  Schoenberg 

For  Leonard  Euler's  bicentenial  of  1983. 

1.  Introduction.  Recently  I  have  written  the  article  [4]  on  this 
subject,  but  the  present  note  can  be  read  independently  of  it.  Its  aim  is  to 
show  briefly  that  the  new  approach  of  §3  leads  to  the  exponential  Euler 
splines  introduced  in  [2].  We  need  a  few  definitions. 

Let  S  =  {s(x)>  (n  >  1)  denote  the  class  of  cardinal  spline  S(x)  of 
n  =  ■ — — — — — 

degree  n.  This  means  that  in  each  unit  interval  ( v,  v  +  1 )  (v  6  Z)S(x)  is 
a  polynomial  of  degree  <  n,  with  the  strong  restriction  that 
S (x )  6  Cn"1(R).  We  also  need  the  class  of  midpoint  cardinal  splines 

*  {s(x ) ;  S(x  +  V2  )  6  S  )•  In  this  case  t*le  junction  points  (or  knots) 
between  the  polynomial  components  of  S(x)  are  at  v  +  V2  (v  6  z)»  Actually 
the  action  will  take  place  within  the  class 


(1.1) 


S  = 
n 


S  if 
n 


S*  if 

n 


A  convenient  basis  for  the  class 
splines  which  are  described  as  follows 


n  is  odd  / 


n  is  even  • 

S  is  furnished  by  the  central  B- 
n 

(see  [3r  Lecture  2,  §1]).  We  define 


(1.2) 


M 1  (X) 


{ 


1  if  J/2  <  x  <  V2 

0  elsewhere  , 


and  convolute  it  with  itself  to  obtain 


(1.3) 


M2(x)  -  M1  *  M ^ (x ) 


-  M  it 

elsewhere 


-1 


<  x  < 

s  s 
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Early  in  this  century  it  was  realized  that  S(x)  »  £  c  M  (x-v)  is  a  unique 

representation  of  every  element  of  S_  *  S,#  This  extends  to  S  as  follows: 

11  n 

We  form  the  convolution 


n+1  factors 


0.4) 


Vi (x) 


M.  *  M. 


M,  (X) 


i-T  mjt:1)**- v*v2n+v2>; 


where  x+  *  max(0,x).  This  an  even  bell-shaped  function  which  is  positive  in 

its  support  (-Votn+I)#  Vb  ( n+1 ) ) •  Moreover  M_.,(x)  6  S_f  and  every 
*  *  n+i  n  — ■  *■ 

S(x)  G  S ^  admits  a  unique  representation  S(x)  ■  ][  cvMn+1  and 

_oo 

conversely  [3,  Lecture  2,  §2]. 


Let  t  be  a  constant  such  that 

(1.5)  t  -  |t|ei0,  -»  <  o  <  w,  t  ?  0,  t  +  1  , 

the  objective  being  to  find  all  solutions  S(x)  G  S  of  the  functional 

n 

equation 

0.6)  S(x+1  )  «  tS(x)f  (x  6  R)  . 

From  [2,  Part  I,  $1]  we  need  the  very  simple  result 
!•  Setting 

m 

(1.7)  G  (X)  -  G  <X,t)  -  l  t^M  (x-V)  , 

n  n  m  n+i 

the  most  general  S(x)  6  S  satisfying  (1.6)  is  of  the  form  S(x)  -  C*G  (x). 
■*  n  ■  mm  ■  —  i  —  ■  n 

From  (3,  Lecture  2,  §(4  and  5]  we  also  need 
II*  If_  t  satisfies  (1.5),  then 

m 

(1.8)  G^tx)  *  E  tVMn+1(x-v)  ^  0  for  all  real  x  • 


This  is  evident  if  t  >  0  because  Gn(x)  >  0  for  all  x.  For  complex¬ 
valued  t  (1.8)  follows  from  the  fact  that  the  curve  of  the  complex  plane 

Z  •  Vx>'  (-•  <  X  <  -) 

spirals  "convexly"  around  the  origin  O  without  ever  reaching  it. 
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2*  Euler ^8 ^genera ting ^function .  Our  actual  discussion  starts  with 
Euler*s  generating  function  [1,  Chap*  VII,  §178] 


(2.1) 


.  4  m  A  <x;t) 

t-1  XZ  V  n  n 

- e  *  )  -■  , —  z 

*  L  n 


t-e  0 

which  defines  the  Appell  sequence  of  monic  polynomials  An(x;t)  *  An(x). 

These  are  the  Euleriah  polynomials*  If  we  differentiate  (2.1)  with  respect 

(v) 

to  x  v  times  (v  >  0)  we  obtain  from  Ay  (x;t)  *  vi  that 


t-1  xz  v 
- -  e  z 

t-eZ 


-  A(V) 

l  -± 


<x»t) 


n*v 

v 


nl 


n  v 
z  =  z  + 


l 

n=v+1 


A^V)(x»t) 
n  n 

ni  Z 


Canceling  the  factor  z  and  setting  successively  x  -  0  and  x  =  1  we 
obtain 


t-l  .  *  »‘v,<o,t)  t_,  , 

- m  I  — - - -  z  ,  -  e  -  1  = 

t-ez  n-v+1  n  t-eZ  n=v+l 


"  An  (1,t)  _n-v 
^  nl  Z 


Multiplying  the  first  equation  by  t,  we  find  that  the  two  left  sides  become 
identically  equal*  Therefore  the  right  sides  are  identical  and  show  that 
(2*2) 


A*W)(1»t)  -  tA*V)(0;t)  (V  «  0,1,...,n-1)  . 

n  n 


This  is  a  characteristic  property  of  the  Eulerian  polynomial  An(x;t). 

The  definition  of  the  exponential  Euler  spline  sn(x)  as  1  defined  it  in 
(2]  is  now  immediate i  We  set 

(2*3)  S  (x )  -  A  (x; t )/A  ( 0? t )  if  0  <  x  <  1  , 

n  n  n  —  * 

and  extend  it*  definition  to  all  real  x  by  the  functional  equation 

(2.4)  Sn(x+1)  -  tsn(x)  (X  9  R)  . 

From  (2.4)  and  the  boundary  property  (2.2)  we  find  that  sn(x)  6  c”"1!*)  and 


therefore 

(2.5) 

Pro*  (2.3)  and  (2.4)  we  obtain  that 

(2.6)  *n(v)  -  tV 


a  (x)  6  S 
n  n 


(v  g  z) 
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(2*7) 


The  exponential  Guler  spline  s  (x)  of  S  is  now  obtained  as 
"  —  1  ■  n  —  n  —  ■  ■  ■  * 

8*(X)  -  8  (X  +  4-)/S  (-i)  (X  «  R)  . 

n  n  2  n  2 

It  satisfies 

(2.8)  s*(x+1)  *  ts*(x)  (x  G  R) 

n  n 

and  therefore 

(2.9)  s*(v)  -  tV  (V  6  2)  . 

n 

*3  ^(Oyt)  f  0  in  (2*3)?  It  is#  for  if  An(0>t)  -  0  we  could  define  a 
cardinal  spline  s(x)  by  setting  s(x)  -  AR(x;t)  in  [0,1)  and  extend  it  to 
satisfy  s(x+1)  *  ts(x).  By  I.  it  follows  that  s(x)  =  C*Gn(x)  (C  f  0), 
and  now  s(0)  -  An(0;t)  -  0  would  contradict  II. 


3*  ^f-ggw^figroacl^to^xgonential  Euler  splines.  As  we  are  going  to 

alternate  between  the  classes  S  and  S*  we  find  it  convenient  to  define 

n  n 


(3.1) 


if  n  is  odd 


S  if  n  is  even 
n 


The  corresponding  splines  we  denote  by 


(3.2) 


s  (x) 
n 


C  8n(x)  if  n  is  odd  # 
l^s^fx)  if  n  is  even  • 


We  now  define  a  new  sequence  of  functions 
(3*3)  Sn(x)  »  Sn(xit)  (n  «  1 # 2# • • • )  # 

by  first  setting 

(3.4)  S1  (x)  -  s^x)  «  s^x)  ■  l  tVM2(x-v)  , 

and  determine  the  functions  (3.3)  by  the  recurrence  relation 

x+l  1 

(3.5)  Sn(x)  -  /  ‘  Sn  (u)du//2  S  (u)du  (n  -  2,3,...)  . 

X“  2  ‘2 
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Our  main  result  is 


Theorem  1.  The  sequence  of  functions  (3*3)  defined  by  (3*4 

identical  with  the  exponential  Euler  splines  of  the  classes  S !  , 

(3.6)  S  (x)  =  s  (x),  (n  *  1 , 2, • • • )  • 

n  n 

Let  us  first  establish  for  the  functions  (1.7)  the  relation 


'  1 

x+  T 

(3.7)  G  (x)  -  /  7  G  (u)du  (n  -  2,3,...)  . 

n  i  n— i 

X-  2 

1 

X+  ~2 

Indeed,  from  (1.4)  we  have  Mn+i<x)  *  M1  *  Mn(x)  *  /  1  Mn(u)du 


1 

rX+  2 

/  M  (u-v)du  -  M  (x-v)  . 

in  n+i 

X-  2 


Now  (1.7)  implies 


1 

x+  — 


1 

x+  — 


/  \  G  (u)du  *  l  tV  /  ?  M  (u-v)du  *  J  tVM  (X-V)  a  G 

1  n-1  ..  'In  “  n+1  n 


X-  2 


X-  2 


proving  (3.7). 


We  observe  next  that  M  (x)  6  S  ;  also  that  by  I*  and 

n+1  n 


(3.8) 

However,  from  (3.7)  we  obtain 


sn(x)  -  Gn(x)/Gn(0)  . 


(3.9) 


V°> 


2 


and  therefore  by  (3.7)  and  (3.3) 


snU) 


X+  2 


f  *  Gn_1(u)du/J  1  Gn_1(u)du 
"  1 


X"  2 


and  (3.5)  is 
hence 


whence 


(u)  , 

II.  we  have 


If  we  define  here  the  terms  of  the  fraction  on  the  right  by  Gn_1 ( 0)  we  find 
from  Gn_.j  (u)/Gn_.j  (0)  •  .^(u)  that 


rx+ 2“ 


( 3e  10)  s  (x)  ■  /  f  S  (u)du/j  s  (u)du  (n  »  2,3,...) 
n  i  n-i  i  n- i 


Since  s«(x)  *  S,(x),  the  relations  (3.10)  clearly  imply  that  S  fx)  ■  s  (x) 
11  «  n 

and  Theorem  1  is  established. 

Remarks.  1.  Notice  that  (3.4)  and  (3.5)  have  produced  only  one  half  of 

the  exponential  Euler  splines,  i.e.  only  those  described  by  (3.2).  If  we 

* 

retain  the  relations  (3.5),  but  start  from  Sj(x)  -  s^(x),  we  would  get  the 

* 

other  half:  s  (x)  ■  s  (x)  if  n  is  odd  and  Sn(x)  •  s_(x)  if  n  is  even. 

n  n  n  n 

2.  If  we  retain  (3.4),  but  modify  (3,5)  to 

Sn(x)  “  ^  Sn-1(u)du/^0  Sn-1(u,du'  (n  ■  2,3,...)  , 

we  would  obtain  that  sn(x)  »  sn(x)  for  n* 

3.  So  fare  we  have  excluded  the  case  when  t  is  real  and  negative. 
Actually,  the  case  t  <  0  leads  to  the  so-called  eigensplines  of  the  classes 

a  n  * 

3^  and  o^,  which  are  fundamental  for  the  problem  of  cardinal  spline 
interpolation  (see  [3,  Lecture  4]).  Exceptional  is  the  case  of  the  Euler 
splines  which  arise  if  t  »  -1 •  In  this  case  our  Theorem  1  again  holds,  as 
shown  in  [4,  §6]. 

4.  In  [2,  §11]  and  again  in  [4, §5]  I  have  recommended  the  use  of  the 
exponential  Euler  splines  for  the  programming  of  the  exponential  function  on  a 
computer.  Especially  for  t  *  2  this  would  produce,  by  a  simple  algorithm, 
very  smooth  and  close  approximations  to  2X. 

5.  We  hope  to  have  shown  in  this  note  the  fundamental  nature  of  Leonard 
Euler's  contribution  to  cardinal  spline  interpolation. 
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